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ABSTRACT 

We use iV-body simulations to model the 12 Gyr evolution of a suite of star clusters 
with identical initial stellar mass functions over a range of initial cluster masses, sizes, 
and orbits. Our models reproduce the distribution of present-day global stellar mass 
functions that is observed in the Milky Way globular cluster population. We find 
that the slope of a star cluster’s stellar mass function is strongly correlated with the 
fraction of mass that the cluster has lost, independent of the cluster’s initial mass, 
and nearly independent of its orbit and initial size. Thus, the mass function - initial 
mass relation can be used to determine a Galactic cluster’s initial total stellar mass, 
if the initial stellar mass function is known. We apply the mass function - initial mass 
relation presented here to determine the initial stellar masses of 33 Galactic globular 
clusters, assuming an universal Kroupa initial mass function. Our study suggests that 
globular clusters had initial masses that were on average a factor of 4.5 times larger 
than their present day mass, with three clusters showing evidence for being 10 times 
more massive at birth. 
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1 INTRODUCTION 


Globular clusters are old, (~ 10 — 12 Gyrs 

dMarin-Franch et al.l 120091 ')') spherical, gravitationally- 
bound configurations of stars with total present-day masses 
in the range ~ 10 4 -10 6 Mq. These clusters are believed to 
have formed within the high pressure environments of the 
primordial Milky W ay, and mark some of the initial sites 
of star formation dKruiissenl I2014T ). Having been found in 
all types of galaxies, we know they must play a pivotal role 
in the galaxy formation process in general. Once a cluster 
forms, stellar evolution and dynamics result i n many of the 
cluster’s initial properties changing over time (iGnedin et al.l 


19991: [F all fc Z hand l200ll: iGieles fe Baumgardtl 20081 : 


Webb et al. 2014a : Brockamp et al.ll2014l . e.g.l. Star cluster 


simulations allow us to trace the evolution of clusters over 
their entire lifetime, and let us compare their present day 
properties to their primordial state. Understanding the 
initial conditions under which a globular cluster population 
forms allows us to explore how galaxies form and further 
constrain galaxy formation models. 

Studies regarding the dynamical evolution of a globular 
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cluster in isolation date back to iHeno 3 (Il96ll~) . After an ini¬ 
tial phase of mass loss via stellar evolution, the dynamical 
evolut ion of star c lusters is driven by two-body relaxation 
fe.g. llTe norJ 1 1 Ufi ll. 1 19731: ISnitzerl 119871 : iHeggie fc Hutl 120031 : 
IGieles. Heggie fc Zhadl201ll '). Relaxation drives mass segre¬ 
gation, causing high-mass stars to sink deeper in the cluster 
potential and low-mass stars to drift outward. In terms of 
structure, the cluster will slowly expand and dissolve, how¬ 
ever the timescale for dissolution in the absence of a tidal 
field is much greater than a Hubble time. The presence of 
a tidal field accelerates the dissolution time of a cluster, as 
two-body interactions only need to push stars beyond the 
tidal boundary as opposed to speeds greater than the clus¬ 
ter’s escape velocity. The tidal boundary will fluctuate over 
the course of a cluster’s orbit if it is non-circular or inclined. 

In a realistic gravitational potential, clusters will also 
be subject to tidal shocks which inject energy into the clus¬ 
ter and temporarily minimize the tidal boundary during 
perigalactic passes, passages through spiral arms and the 
Galactic disk, a nd encounters with gian t molecular clouds 
(GMCs). In fact. lLamers fc Gielesl l2006i l finds that encoun¬ 
ters with GMCs is the dominant mechanism behind evapora¬ 
tion for clusters orbiting in the solar neighborhood. Simula¬ 
tions bv lKruiissen et akl (l201lh of clusters with and without 
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tidal shocks illustrate that tidal shocks will lead to between 
80 — 85% additional dissolution events depending on the 
galaxy model. While tidal shocks experienced by clusters 
within the high-pressure environment in which they form 
would lead to very short dissolution times, the distribution 
of clusters in the Milky Way indicate that clusters survive by 
migrating outwards into the halo. The migration can likely 
be attributed to galaxies forming and growing via the h icrar- 
chical merger of smaller galaxies (ISpringel et al.ll2005l). with 
merg er events push ing clusters out to wide r orbits llKruiissenl 
120141) . However as iKruiissen et al.l (I2012T) points out, tidal 
shocks via merger events also have the ability to significantly 
disrupt clusters, sometimes to the point of complete dissolu¬ 
tion. If a cluster survives to reach its present day orbit in the 
halo, the tidal shocks experienced by the cluster will be much 
weaker than when it was in its formation environment, with 
the dissolution times of surviving clusters becoming greater 
than a Hubble time. 

When a cluster is stripped of stars, regardless of whether 
the mass loss event is due to two-body interactions, tidal 
stripping, or tidal shocks, it will primarily be stars orbiting 
in the outskirts that are removed. Hence over time, since 
two-body relaxation leads to mass segregation, the proba¬ 
bility that a given star will escape from a cluster into the 
tidal field of its host galaxy over the course of a relaxation 
time increases as a function decreasing mass. In turn, the 
mean stellar mass in clusters slowly increases over time, and 
with a rate that decreases with increasing total cluster mass 
since the time-scale for two-body relaxation scales as N 1,/2 , 
where N is the number of objects in the cluster. 

IVesperini fe Heggiel (119978 1 conducted one of the first 
studies that explored the effects of different evolutionary 
processes on a cluster’s stellar mass function. The authors 
considered the effects of relaxation, stellar evolution, disc 
shocking a nd a tidal field on globular c luster evolution. A key 
result from IVesperini fc Heggiei (Il997l f is that the difference 
between the initial and present day mass function of a cluster 
is primarily dependent on the fraction of mass lost, with 
only minor scatter due to initial mass, initial concentration, 
and Galactocentric dista nce (assuming circ ular orbits). This 
findin g was confirmed bv lKrui isscn (20093) and lTrenti et al.1 
(2010), however the former noted that the retention fraction 
of black holes could alter the evolution of the mass function. 
More specifically, a cluster hosting either a central massive 
black hole or multiple stellar-mass black holes will have a 
higher escape rate of high-mass stars due to their interaction 
with black holes in the cluster core. The mass function will 
in turn evolve at a much slower rate as lo w-mass stars are no 
longe r the only stars esca ping the cluster l|Lutzgendorf et ahl 
|20l.‘j) . [hC r 11 i isseiJ (l2009bl ) also finds that the evolution of the 
mass function can be dependent on the dissolution time of 
the cluster, as the mass function of young clusters with short 
dissolution times evolves slowly if most of the mass loss is 
occurring before high mass stars reach the end of their life. 

iWebb et al.l ( 2014al ) builds on previous work and ex¬ 
plored the effects of orbital eccentricity on the evolution 
of the _stellar mass function slope a. The main result of 
iWebb et al.l (i2014a ) is that for a given perigalactic distance 
R p , the dynamical evolution of the mass function slows with 
increasing orbital eccentricity due to the weaker mean tidal 
held it experiences. Thus the slope of the mass function 
a increases (i.e. becomes less negative) over time, with a 


rate that increases as the mean strength of the tidal held 
increases. When comparing model clusters with the same 
fraction of initial mass lost, the evolution in a was almost 
orbitally-independent. For a given fraction of mass lost there 
was only a minor spread in a due to cluster orbit as stronger 
tidal helds can remove higher mass stars from the cluster 
(i.e. the dependence of the probability of escape on stellar 
mass becomes slightly weaker in a stronger tidal held). This 
result further strengthens the possible relationship between 
the mass function of a globular cluster at a given time and 
the fraction of its initial mass t hat has been lost. However, 
we noted in i Webb et al.l (l2014al ) that differences due to the 
combined effects of orbital eccentricity, orbital inclination, 
initial size, and initial mass of a globular cluster may in some 
cases slightly alter the relation. 

The present paper is aimed towards quantifying how 
these combined effects will alter the evolution of a as a 
function of fraction of initial mass lost. lLamers et al.l (120131 ) 
notes that comparing clusters as a function of initial mass 
lost as opposed to time ensures that clusters are being com¬ 
pared on the same evolutionary timescale. If a globular clus¬ 
ter’s stellar mass function can be shown to be primarily de¬ 
pendent on the fraction of mass lost, it will also represent 
the clusters’s dynamical age. Hence the initial mass of a 
cluster can be determined assuming that all globular clus¬ 
ters formed with the same stellar initial mass function (IMF) 
and either: 

• A) tidal shocks experienced by clusters before they mi¬ 
grate to their present day orbit remove stars over the en¬ 
tire mass spectrum, such that the mass function does not 
evolve. Hence all clusters will still have the same IMF once 
they reach the halo and their mass function will reflect their 
dynamical age since escaping from their formation environ¬ 
ment. If this is the case, the initial cluster mass which cor¬ 
responds to the cluster’s present day mass function will rep¬ 
resent its mass at the time of migration into the halo (i.e. 
its current orbit). Or, 

• B) tidal shocks experienced by clusters before they mi¬ 
grate to their present day orbit remove stars in a similar 
manner as tidal evaporation and disk shocking, just over a 
shorter time scale. Hence the evolution of the mass function 
will remain unchanged with respect to the fraction of mass 
lost, and the mass function will reflect a cluster’s dynamical 
age since formation. In this scenario, the initial cluster mass 
calculated from the cluster’s present day mass function will 
represent its true initial mass at formation. 

However if tidal shocks that occur while clusters are still 
within their formation environment remove stars such that 
the evolution of the mass function is affected in a manner 
different from the above picture, then the connection be¬ 
tween a cluster’s stellar mass function and its dynamical age 
weakens. We will henceforth refer to this kind of mass loss as 
’’non-standard”, since the underlying physical mechanisms 
are unknown. 

The assumption that all clusters form with the same 
IMF is consistent with the observed present-day stellar 
mass functi ons, once thei r dyn amical evolution has been ac¬ 
counted for llLeigh et al. Il20l3) . Hence the possibility of sim¬ 
ply using the mass function of present day globular clusters 
to predict their initial masses becomes an exciting possibil¬ 
ity. How exactly globular cluster populations evolve from 
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their initial cluster mass function (GCIMF) to the present- 
day observed globular cluster mass function (GCMF) is still 
a topic of much debate. 


The GCMF has been found to be nearly universal, 
a gaussi an centred around a m ean mass of approximately 
10 5 Mq llBrodie fc Straderll200fil l. Whatever its initial form, 
the GCIMF will have undergone significant evolution due 
to clusters being sub jecte d to tidal shocks within the for¬ 
mation environment llKruiissenll2014l ) , internal two-body re¬ 
laxation, dyna mical friction, disk shocks and bulge shocks 
d Gnedin et al.1119991 ). It is not clear whether the present 
day mass function contains any traces of the GCIMF. The 
GCIMF is often assumed to have a power-law form, as 
this would would match the mass functions of other as¬ 
tronomical objects like molecular clouds and young clus- 
tersjRdM^^Zhang 2001). A similar conclusion was reached 
by iMcLauehlin fc Falll ( 20081 ) , who finds that the Galac¬ 
tic GCMF depends on cluster half-mass density, which 
they argue is a signature of a mass function that ini¬ 
tially rose towards low masses but has eroded over time 
due to internal two-body relaxation driving the dissolu¬ 
tion of its (preferentially low-mass) star clusters. How¬ 
ever both these early studies make the simplifying as¬ 
sumption tha t cluster mass-loss rat e s are environment ally- 
independent. iGieles fc Baumgardtl d2008l) argued that a 
power law GCIMF will not evolve to become the present-day 
GCMF due to internal two-body relaxation alone, whereas 
other analytic prescriptions and globular cluster models 
have been able to reproduce the present-day GCMF by 
including not only the effects of evaporation by two-body 
relaxation, but also gravitational shocks due to a tidal 
field, stellar evolution , dynamical friction, gas expulsion 
and radial aniso tropy ( Fall fc ZhandlioQll : IVesnerinill200ll : 
iBrockamp et al.ll2014l ). Other studies have also been able 
to reproduce the present-day GC MF beginning with a log¬ 
normal GCIMF (|Vesperini 111998 )1. All of these studies are 
based on treating the present-day GCMF and the GCIMF 
as probability distribution functions which can be altered 
due to several processes. A method for determining the ini¬ 
tial mass of the remaining globular clusters in the Milky Way 
would provide robust constraints on the possible formation 
of globular cluster systems and galaxies themselves, while 
also providing insight into whether we understand how the 
various mechanisms discussed above alter the GCIMF over 
time. 


In t his stud y, we combine the results presented in 
IWebb et al.1 (l2014al ) with additional models that range in 
initial mass, initial size and orbit in order to develop a 
method of predicting the initial total stellar mass of a glob¬ 
ular cluster based on its present-day stellar mass function. 
The effects of two-body relaxation are quantified by consid¬ 
ering a range of cluster masses and sizes. By placing our 
model clusters on a range of orbits we also consider the ef¬ 
fects of tidal heating, tidal shocks, and the initial degree 
of tidal-filling. We also allow stars in our model cluster to 
evolve, hence stellar evolution is also being considered. In 
Section [2] we introduce the models used in this study and in 
Section [3] we compare the value of a for each model cluster 
at different time steps to the fraction of mass each cluster 
has lost. In Section [4] we will quantify the relationship be¬ 
tween a and the fraction of mass lost by a cluster and use 


it to predict the initial masses of select Galactic globular 
clusters. We summarize our findings in Section [5] 


2 N-BODY MODELS 


The direct Wbody code NBODY6 dAarseth et al.l Il974l : 
lAarseth~ll20 03) is used to model star clusters o ver 12 Gyr 
of evolution. We use a P lummer density profile (PI u mmer I 
1 191 ll : lAarseth et al.lll974l) to generate the three dimensional 
position and velocity of each star within the cluster out to a 
cut-off of ~ 10 r m , where r m is the system half-mass radius. 
Since we are primarily concerned with the evolution of the 
stellar mass function, the stars within each model cl uster 
are as signed masses based on a iKroupa. Tout, fc Gilmo~ 
(11993h IMF (i.e. the IMF found by ILeigh et al. i ( 2012) 
that best reproduces the observed cluster-to-cluster varia¬ 
tions in the present-day stellar mass functions in a large 
samp le of Galactic GCs) betwe en 0.1 and 30 Mq. The 
iKroupa, Tout, fc Gilmor e (jl993l ) IMF has the functional 
form: 


dN 

dm 


a 


= m 


(i) 


where a equals -2.7 for m > 1 Mq, -2.2 for 0.5 ^ m ^ 
IMq, and -1.3 for 0.08 < m ^ 0.5Mq. With respect to 
Equation [lj the comm only cited Salpeter IMF has an a of 
-2.35 dSaloeter II 19551 '). 

Model clusters all have metallicities of Z = 0.001 and bi¬ 
nary fractions of 4%. Binary stars are assigned total masses 
based on two stars drawn from the IMF, while the mass 
ratio between the two stars are selected from a uniform dis- 
tributionAThe initial dis tribu tion of binary periods matches 
that of lDuquennov fc Mavod dl99ll) while the initial distri¬ 
bution of orbital eccentricities follow a thermal distribution 
dHeggie ill975l) . Details regarding the stellar and binary evo¬ 
lution algorithms us ed in our simul ations can be found in 
iHurlevI (|2008al ) and iHurlevI (|2008bl ). In order to compare 
the results of our models to Galactic globular clusters, our 
model clusters orbit within a Milky Way-like potential com- 
prised of a 1.5 x lO^Mp, point-mass bulge, a 5 x 10 10 Mq 
iMivamoto fc Nagail d 19751) disk (with a = 4.5 kpc an d b = 
0.5 kpc), and a logarithmic halo potential dXue et al.ll2008l ). 
The logarithmic halo potential is scaled such that the cir¬ 
cular velocity at a galactocentric distance of 8.5 kpc is 220 
km/s. The majority of mo dels discussed in th i s study were 
initial ly p resented in either I Webb et aid d2013l ) , ILeigh et al. I 
d2013l) , or I Webb et~aD d2014bl) . Additional models only dif¬ 
fer in initial mass and/or size. We refer the reader to these 
previous studies for a more detailed description of our sim¬ 
ulations. 

Our initial suite of simulated clusters had initial masses 
of either 3 x 10 4 Mq or 6 x 10 4 Mq and initial half-light radii 
rh of either 2 pc or 6 pc. The models are taken from pre¬ 
vious works studying the effects of initial size and orb i t on 
the dynamical evolution of star clust ers (IWebb et al.ll2013l ; 
ILeigh et al. Il201~3l : IWebb et al.ll2014al ). To further study the 
effects of initial mass on cluster evolution, select clusters 
were re-simulated but with initial masses of 8 x 10 4 Mq and 
1.1 x 10 5 Mq. For a given initial mass and size, we model 
clusters with orbital eccentricities of 0, 0.25, 0.5, 0.75, and 
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Figure 1. Illustration of star cluster orbits that are used in 
this study. Each initial mass and size combination was modelled 
with a circular orbit at 6 kpc (solid black line). Additional orbits 
that correspond to various initial mass and size combinations are 
marked in the legend. For inclined orbits, star clusters were mod¬ 
elled with orbital inclinations of 22° and 44°. 

0.9, where eccentricity is defined as e = ■ Ra and R p 

are the apogalactic and perigalactic distances of the orbit, 
respectively. All models with eccentric orbits have a R p of 6 
kpc. For comparison purposes, a cluster was also modelled 
with a circular orbit at the R a of each eccentric cluster. 
From a previous study focusin g on the effects of o rbital in¬ 
clination on cluster evolution rtWebb et al.| [ 2014bfl . we also 
include model clusters with orbits having inclinations of 22° 
and 44° relative to the plane of the Galaxy. These additional 
models also include clusters with inclined and eccentric or¬ 
bits. Since our previous model clusters with inclined orbits 
all had an initial i'h of 6 pc, we have re-simulated these clus¬ 
ters with an initial th of 2 pc to explore the effects of initial 
size on cluster evolution. The initial model parameters of 
our entire suite of simulations are summarized in Table [T] 
The orbits covered by our suite of simulations are illustrated 
in Figure [l] 0 


3 THE STELLAR MASS FUNCTION - INITIAL 
MASS RELATION 

To find the value of a in Equation |T] for each model clus¬ 
ter at a given time step, we first mark all main sequence 
stars within the tidal radius of the cluster that have masses 
between O.15M0 and 0.8 Mq. The latter mass corresponds 
approximately to the main sequence turn off at 12 Gyr. Ex¬ 
cluding main sequence stars beyond the main sequence turn 
off has a negligible effect on the slope of the mass function, 

1 It should be noted that all model clusters have orbital velocities 
greater than 88 mph. 


Table 1. Model Input Parameters 


Mi 

fh,i 

Rp 

e i 

Reference 

Mq 

pc 

kpc 

degrees 



3 X 10 4 2 6 0 0 Webb et al. ('2013') 



6 

0.25 

0 



10 

0 

0 



6 

0.5 

0 



18 

0 

0 



6 

0.75 

0 



43 

0 

0 



6 

0.9 

0 



104 

0 

0 


3 x 10 4 6 

6 

0 

0 

Webb et al. (20131 


6 

0.25 

0 



10 

0 

0 



6 

0.5 

0 



18 

0 

0 



6 

0.75 

0 



43 

0 

0 



6 

0.9 

0 



104 

0 

0 


6 x 10 4 6 

6 

0 

0 

Leigh et al. (2013) 


6 

0.5 

0 



18 

0 

0 



6 

0.9 

0 



104 

0 

0 


2 

6 

0 

22 

This publication 


6 

0.5 

22 



18 

0 

22 



6 

0.9 

22 



104 

0 

22 



6 

0 

44 



6 

0.5 

44 



18 

0 

44 



6 

0.9 

44 



104 

0 

44 


6 

6 

0 

22 

Webb et al. (2014b) 


6 

0.5 

22 



18 

0 

22 



6 

0.9 

22 



104 

0 

22 



6 

0 

44 



6 

0.5 

44 



18 

0 

44 



6 

0.9 

44 



104 

0 

44 


8 x 10 4 6 

6 

0 

0 

This publication 


6 

0.9 

0 


1.1 x 10 5 6 

6 

0 

0 

This publication 


6 

0.5 

0 



since these phases of stellar evolution represent a very nar¬ 
row range in stellar mass. This exclusion also ensures that we 
are not including main sequence stars that have formed via 
collisions or mergers. In order for our models to be compared 
to observations, we assume binary stars are unresolved and 
treat them as single stars. The total mass of the binary is 
calculated from the total luminosity and assuming L~M5. 
We then set a equal to the slope of a plot o flog(^-) versus 
log(m). It should be noted that since the mass function is 
best represented by a two-part power law, we take a to be 
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Figure 2. Slope of the stellar mass function at 6, 10, and 12 Gyr 
versus the fraction of mass lost for each model cluster. Points 
are colour coded based on initial mass while different symbols 
represent different initial sizes as indicated. Exponential curves 
of best fit are matched to all data points (solid line), clusters 
with initial sizes of 6 pc (dotted line) and 2 pc (dashed line). 


the slope of the high mass end of the mass function, beyond 
O.15M 0 . 

To get a sense of how the mass function of each model 
cluster evolves, we plot a for each model at 6, 10, and 12 
Gyr as a function of fraction of mass lost in Figure [2] It 
should be noted that some of the lower initial mass models 
reach dissolution before 12 Gyr, and therefore only have two 
data points in the figure. 

As expected, the general trend is for a to increase as a 
cluster loses mass due to either two-body interactions eject¬ 
ing stars from the cluster or tidal stripping and bulge and 
disc shocks removing primarily low mass sta rs that have 
been mass segregated to the outer regions (e.g. iLeigh et al. I 
120121 ). In order to track the evolution of a as a function of 
mass lost, we fit our data points with an exponential curve 
of the form: 

R M 

a = A x e ^ + C (2) 

In the absence of stellar evolution and binary stars, the 
offset C would simply be the initial slope of the mass func¬ 
tion on. Since we include the effects of stellar evolution and 
binaries here, the offset C must compensate for a initially 
decreasing due to the cluster losing mass via st ella r evo lution 
and the disruption of binary stars rtWebb et al.ll2014al ). Dur¬ 
ing this phase, which lasts until clusters have lost approx¬ 
imately 40% of their initial mass, clusters will also expand 
and lose some mass in the form of stars. For our models, C 
is set to be just larger than the most negative value of a. 
Since we do not model the initial decrease in a, our relation 
is only applicable to star clusters that have lost greater than 
40% of their initial mass. This criterion is met by Galactic 


globular clusters which will have lost a significant amount 
of their initial mass from stellar evolution alone. 

Initially taking into consideration all models, we find 
that A = —2.7±0.1 and B = —6.4±0.1, which is plotted as 
a solid curve in Figure[5] Here we have established an evolu¬ 
tionary track for a that is solely d ependent on the a mount 
of mass a cluster has lost. However. iKruiissenl ll2009bl ) notes 
that the evolution in a will be different for clusters with 
short d issolution times. Compar isons of Equation [2] to the 
work of lLutzgendorf et al.l d2013l ) suggests that the existence 
of a massive black hole (greater than 1% of the total clus¬ 
ter mass) at the centre of a cluster may alter our predic¬ 
tion of jf- by up to 20%. Observations of such massive 
black holes in glo b ular clusters have yet to be confirmed. 
iLutzgendorf et all 1 2013 ) also finds that the black hole re¬ 
tention fraction can alter the evolution of a, however the 
added uncertainty in jj. of ~ 10% is comparable to the 
scatter about Equation [2] as discussed below. 

The scatter about the relation given by Equation [2] is 
only a = 0.09, and is primarily due to cluster orbit and ini¬ 
tial cluster size, while initial cluster mass appears to have 
no affect on the evolution of a. We recover the actual initial 
mass of each of our models to within 7% when usin g Equa¬ 
tion [2] As previously found in IWebb et alj J2014al) . cluster 
orbit can slightly alter the evolution of a since stronger tidal 
fields can remove higher mass stars than weaker tidal fields. 
When considering clusters of different initial size, a larger 
amount of scatter is introduced to the mass function - ini¬ 
tial mass relation. An under-filling cluster will evolve as if it 
is in isolation, undergoing a certain degree of two-body re¬ 
laxation and mass segregation until it expands to the point 
of becoming tidally affected (fG ieles. H eggie fo Z hao 2011; 
lAlexander fc Gielesl 120131 : 1 Webb et alJ 12013 ) . This E quali¬ 
tatively equivalent to assuming some degree of primordial 
mass segregation in a tidally-filling cluster. When an ini¬ 
tially under-filling cluster begins losing stars it is preferen¬ 
tially losing low mass stars that have been segregated to the 
outer regions, which is reflected in a faster increase in a as 
a function of the fraction of mass lost. Since a tidally filling 
cluster is subject to tidal stripping almost immediately, a 
mix of high and low mass stars is removed before segrega¬ 
tion becomes important, resulting in a shallower increase in 
a initially. 

If we choose to fit the extended (circles - rh,i = 6 pc) and 
compact (squares - rn,i = 2 pc) clusters separately, we End 
that for extended clusters A = — 2.4±0.1 and B = —6.6±0.1 
(plotted as a dashed line). With the exception of the clus¬ 
ters orbiting beyond 18 kpc, all extended models are tidally 
filling at R p and for a higher fraction of their orbit if their 
eccentricity is low. Hence these clusters can lose stars over a 
higher mass range initially, slowing the evolution of a com¬ 
pared to the overall population. For compact clusters we find 
A = —2.8 ± 0.1 and B = —5.8 ± 0.1 (plotted as a dotted 
line). All of the compact clusters are initially tidally under¬ 
filling, such that they have time to undergo some degree 
of mass segregation before tidal stripping becomes a major 
source of mass loss. Hence the evolution of a is faster than 
the overall population since primarily low-mass stars are be¬ 
ing removed. It should be noted that the uncertainty in the 
extended and compact best fit relations are 5% and 2.5% 
respectively, indicating that knowing a cluster’s initial size 
will only marginally improve the determination of its initial 
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Table 2. Model Input Parameters 


Dataset 

A 

B 

C 

0.15 - 0.8 Mq 
All 

2.8 

-6.4 

-1.5 

Compact 

2.8 

-5.8 

-1.5 

Extended 

2.4 

-6.6 

-1.5 

0.1 - 0.5M© 

All 

2.1 

-5.7 

-1.4 

Compact 

2.2 

-5.3 

-1.4 

Extended 

1.9 

-5.7 

-1.4 

0.3 - 0.8M© 
All 

3.8 

-6.7 

-1.8 

Compact 

3.9 

-6.0 

-1.8 

Extended 

3.2 

-6.8 

-1.8 

0.5 - 0.8M© 
All 

4.6 

-7.2 

-2.2 

Compact 

4.6 

-6.4 

-2.2 

Extended 

4.0 

-7.4 

-2.2 


mass via the mass function - initial mass relation presented 
here. However as we will discuss in Section[4] the uncertainty 
due to the omission of additional mass loss mechanisms from 
our models is expected to dominate compared to the uncer¬ 
tainty in initial cluster size. 

We draw the reader’s attention to the specific set of ex¬ 
tended data points that are significantly removed from the 
dashed line of best fit in Figure [2] These open circles rep¬ 
resent tidally-Hlling clusters with circular orbits at 6 kpc 
b ut wit h or bital in clinations of 22° and 44°. As found in 
[Webb et al.1 (l2014blh clusters with inclined orbits at low R gc 
are subject to significant amounts of tidal heating and disk 
shocking which result in significant stellar mass loss. Due 
to the strength of the tidal field and the fact that these 
clusters are not initially mass segregated, higher mass stars 
can be removed from the cluster than if its orbit was not 
inclined, slowing the evolution of a. Since these clusters de¬ 
viate slightly from our derived mass function - initial mass 
relation, we suggest that if the orbit of a cluster is known to 
be inclined and at low Galactocentric radius, the line of best 
fit found by considering only our initially extended models 
will yield the initial cluster mass to within a few percent. 

The best fit values to both the entire dataset as well as 
the extended and compact clusters are listed in Table [2] In 
order to apply our relation to a wider range of observational 
datasets, we also include in Table [2]the best fits values of A, 
B, and C when only stars within 0.1 — 0.5Mq, 0.3—0.8 Mg, or 
0.5 — 0.8 Mq are used to find a. Using the lower mass range 
our relation recovers the initial mass of each model cluster 
to within 8.6% and using the high mass range recovers the 
initial mass of each model cluster to within 7%. 


4 DISCUSSION 

The near universality in the evolution of a suggests that as a 
cluster loses mass, the same distribution of stellar masses are 
either tidally stripped or removed via two-body interactions 
from the cluster no matter when the escape occurs. This 
statement is however only applicable to star clusters with 
dissolution timescales greater than the time scale for stel¬ 


lar evo lution, which is true for globular clusters (iKruiissenl 
l2009bl 'l. We do find that the strength of the tidal field and 
amount of mass segregation a cluster undergoes before be¬ 
coming tidally affected can influence the stellar mass func¬ 
tion evolution, however the effects are minimal. Therefore, 
it is possible to use the mass function - initial mass relation 
presented here to determine the initial total stellar mass of 
any Galactic globular cluster. 

Our definition of initial total stellar mass is best de¬ 
scribed as the mass that a globular cluster has when its stel¬ 
lar mass function begins to deviate from its initial value. The 
calculation of the initial mass, since it stems from the sim¬ 
ulations presented above, is based on two key assumptions. 
The first assumption is that all clusters form with the same 
IMF. While the evolution of the differential mass function 
has b een shown to be independent of the IMF (iLamers et ali 
l2013h . an IMF would still have to assumed in order to use 
a cluster’s differential mass function to calculate its initial 
mass. Therefore, to test the dependence of our mass func¬ 
tion - initial mass relation on this assumption we present 
additional simulations of our model cluster with a circular 
orbit at 6 kpc, initial mass of 6 x 10 4 M©, and initial half 
light radius of 6 pc with a range of IMFs. More specifically 
we model clusters with pow er law IMFs tha t range in ca, 
with cti = —2.35 equalling a lSalpeterl d 1955l l IMF. To best 
compare models with different IMFs, we plot the evolution 
of Aa = a;; — a for each cluster in Figure[3]with a calculated 
in three different mass ranges. For the lower mass range, we 
also co mpare our results to previous studies by IKruiissenl 
ll2009bl j (K09) and I Vesperini fe Heggiel (Il997l j (VH97). 

Figure [3] indicates that our mass function - initial mass 
relation is not independent of the stellar IMF. Clusters with 
top heavy IMFs are examples of how the evoluti on of a is 
differe nt for clusters with short dissolution times (|Kruiissenl 
l2009bl h as a significant amount of mass is lost quickly due 
to stellar evolution while the low mass end of the mass func¬ 
tion remains unchanged. The early mass loss also leads to 
significant cluster expansion and the short cluster dissolu¬ 
tion time. Hence a stays relatively constant until the cluster 
is close to completely evap orating. 

For clusters with a ISalpeteiH (Il955l j IMF or power 
law IMF with eg = —2.5, which are similar to 

iKroupa, Tout, fe Gilmorel |l99i|) over the mass range in 
which a is calculated, Equation[2]is able to recover the initial 
mass of a cluster to within 25% and 8% respectively when 
stars between 0.15 — 0.8M© are used. Using the lower mass 
range provides initial mass estimates comparable to when 
the full mass range is used, however the higher mass range 
yields errors approximately 10% higher due to its sensitivity 
to stellar evolution. If a cluster has a bottom heavy IMF, 
very little mass loss will occur due to stellar evolution such 
that a begins to decrease earlier than predicted by Equation 

m 

Compariso n of our mass funct ion - initial mass relation 
to the works of lKruiisse~nl ll200 9bj and lYesperini fe Heggiel 
(Il997l) also provides insight into its dependence on the IMF, 
stellar evolution, and the initia l binary fraction. Th e discr ep- 
ancy between our models and lYesperini fe Heggiel (|l997l 'l is 
a bit surprising, as their models have a power law IMF with 
cti = —2.5 and they take into consid erate the effects of 
stellar evolution. Hence the results of IVesperini fe Heggiel 
d 19971) should be directly comparable to our models in Fig- 
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Figure 3. Change in the slope of the stellar mass function from 
its initial value versus the fraction of mass lost for model clus- 
ters with different IMFs. Functional forms of the IMFinclude 
iKroupa, Tout, Sz Gilmorel \ 19931) (solid black line), ISalpeterl 
(|1955l ) (blue), and power laws with each Qj provided in the legend. 
At each time step a is calculated using stars within 0.1 — 0.5 Mq 
(left panel), 0.15 — 0.8 Mq (centre panel) and 0.3 — 0.8 Mq (right 
panel). In the le ft panel, th e evolution of a as determ ined by 
IKruiissenl (|2009bh (K09) and IVesperini & Heggiel <|l997t) (VH97) 
are also included. 


ure [3] The key differences between the two studies are that 
our maximum stellar ma ss (50 Mq) is much higher than 
IVesperini fc Heggii (| 19971 ) (15Afg) and our study also in¬ 
cludes an initial binary fraction of 4%. A high maximum 
stellar mass will result in significant early mass loss before a 
evolves from its initial value. And as previously mentioned, 
binary stars also serve to delay the evolution of a as bi¬ 
nary stars are disrupted and single stars population the low 
mass end of the ma s s func tion. The decrease in a found by 
IVesperini fe Heggiel (119971 1 is similar to ours if adjusted for 
the different fractions of mass lost before a begins to deviate 
from its initial value. 

Theor etical m odels for the evolution of Aa given by 
iKruiissenl d2009bl f on the other hand are in strong agree¬ 
ment with our own for -AT ^ 0.5 (see Figure [3]). The main 
difference is only in the early evolution of Aa, which as pre¬ 
viously stated is sensitive to the maximum stellar mass and 
initial b inary fraction. Hence both this study and iKruiissenl 
d2009bl l offer methods for calculating the initial mass of a 
globular cluster based on the change in the slope of the mass 
function from its initial value, assuming -AT ^ 0.5. Further 
simulations are required to determine how Equation [2] de¬ 
pends on the maximum stellar mass and initial binary frac¬ 
tion. 

The second assumption associated with our calculation 
of initial mass stems from the fact that we assume clus¬ 
ters have always had their present day orbit. We have not 
modelled stellar clusters within their formation environment 
(where tidal shocks from GMCs have been shown to be a 


major source of mass loss) or as they migrate from their for¬ 
mation environment to their current orbit. How this early 
stage of globular cluster evolution affects the evolution of 
the stellar mass function is still not well understood, and 
may influence the calculations presented here. However, as 
discussed in Section[T]there remain two possibilities in which 
the effects of early tidal shocks may have a negligible impact 
on our results. First, if these early tidal shocks remove stars 
over the entire stellar mass range, then clusters will more 
or less retain their primordial IMF when they migrate out¬ 
wards to the halo. Hence the initial mass presented here will 
simply reflect the cluster’s mass once it reaches the halo. If 
these early tidal shocks acting on globular clusters soon after 
they form do not remove mass from the cluster in some non¬ 
standard way, our calculated initial stellar mass will then 
effectively represent the cluster’s actual initial stellar mass 
at birth. More specifically, as long as tidal shocks are re¬ 
moving the outermost stars in the cluster the same as tidal 
evaporation, albeit at an accelerated rate, the evolution in 
a as a function of the remaining mass fraction should re¬ 
main the same. If anything, since early tidal shocks are es¬ 
sentially stronger and more frequently occurring versions of 
disk shocks, the evolution of a is likely to be faster than 
presented here. In this case, the initial masses calculated 
with Equation [2] are only lower limits. Detailed simulations 
of clusters in this early phase over a range of initial condi¬ 
tions are required before we can quantify its effect on the 
mass function - initial mass relation presented here. 


4.1 Application to Milky Way Globular Clusters 

Our study on the evolution of a cluster’s stellar mass func¬ 
tion can be applied to observations of globular clusters us¬ 
ing the observed present-day global mass function slope, or 
even the mass function slope within a narrow radial bin 
at rh, since this has b een shown to b e comparable to th e 
global mass function (iD e Marchi, Paresce fe Pulone) |2007|: 
DeMarchi, Paresce fe Porteeies Zwartj _ 2010|; IPaust et all 
2010i). From the datasets of IPaust et ali ( 2010i l and 
De Marchi, Paresce fe Porteeies Zwarfl (|2010|), comb ined 


with the Harris GC catalog ( Harrisl 1996. 2010 updatel) . we 

have the present-day total stellar mass (via the integrated 
V-band magnitude of the clusters) and values of a for 33 
Galactic globular cl usters. For clusters in both datasets, we 
take a values from [Paust et all l(2O10l ). Given Equation [5] 
we estimate the initial total stellar mass of each of these 
clusters and plot their initial and present day total stellar 
masses in Figure 0] Since each a was measured over dif¬ 
ferent mass ranges within 0.1 — 0.8 Mq, we use the most 
appropriate mass range from Table [2] when applying Equa¬ 
tion 0 to each cluster. The red dotted lines in Figure [I] 
mark contours of constant a while open circles represent 
clusters for which we could only calculate an upper limit 
for their initial mass. Table [7] contains the present day pa¬ 
rameters, initial masses and orbital parameters (when pos¬ 
sible) of each cluster in our datas et. Or bital parameters 
are taken from Dmescu_e^_al. ICasetti-Dinescu et al.l 

d2007f l and ICasetti-Dinescu et all (|2013i ). We stress that it 
is only possible to constrain the total initial stellar masses 
in our sample of globular clusters. Assuming some star for¬ 
mation efficiency < 100%, the total initial cluster masses 
(including both gas and stars) will have been even higher. 
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Log (Mj/M 0 ) 

Figure 4. Present day globular cluster mass as a function if initial 
mass calculated via Equation [2] Open circles mark clusters for 
which only an upper limit could be placed on their initial mass. 
Red dotted lines represent lines of constant a equal to -2.0 (lower 
red line), -1.5, -1.0, -0.5, 0.0, and 0.5 (upper red line). 

It should be noted that only upper limits could be 
placed on the initial masses of NGC 6304 and NGC 7078 be¬ 
cause their present day a is less than ct; in the 0.3 — 0.8Mq 
mass range. The simplest explanation for this discrepancy is 
th at the 0.3— 0 .8 Mp, mass r ange e ncompasses the bend in the 
iKroupa, Tout, fe Gilmore! d 19931 ) IMF at 0.5Mg, which can 
strongly affect measurements of a depending on the cluster’s 
dynamical age. To place an upper limit on each cluster’s ini¬ 
tial mass, we use Equation[2]with the 0.5—0.8Mq mass range 
since at is less than each cluster’s present day a. Followup 
measurements of a between 0.1 —0.5Mq or 0.5 —0.8Mq will 
allow for a proper estimation of the initial masses of NGC 
6304 and NGC 7078. 

Figured] suggests that the clusters in our data set were 
on average a factor of 4.5 times more massive than they 
are today, with some clusters having initial masses greater 
than a factor of 10 times (e.g. NGC 2298 and NGC 6838) 
or even 20 times (e.g. NGC 6712) larger than their present 
day values. Based on our limited data set, we estimate that 
the initial globular cluster mass function was centred around 
LOG(jj^) = 6.0. However, it should be stressed that this is 
the initial globular cluster mass function of clusters that sur¬ 
vived to the present day only, and does not include clusters 
that have dissolved du e to relaxation, dyn amical friction, or 
disk and bulge shocks ([Gnedin et a.1.11 19991 ). The inclusion of 
these clusters will change the initial glob ular cluster mas s 
function substantially dKruiissen fe Porteeies Zwartll2009d) . 
and should ultimately contribute to lowering our calculated 
mean mass. 

Given the mass function of a globular cluster, our pa¬ 
rameterization (Equation [7]) provides an estimate of ini¬ 
tial cluster mass as a function of its present day mass 
function slope. However, the implied mass loss is on av- 


Table 3. Initial Masses of Galactic Globular Clusters 


ID 

M 

a 

Mi 

(.logjfe) 

Rp 

(kpc) 

Ra 

(kpc) 

Pal 5 

4.30 

-0.4 

5.13 



NGC 104 

6.00 

-0.84 

6.77 

5.2 

7.3 

NGC 288 

4.93 

-0.83 

5.62 

1.7 

11.2 

NGC 362 

5.60 

-1.69 

6.12 

0.8 

10.6 

NGC 1261 

5.35 

-0.59 

6.12 



NGC 1851 

5.56 

-0.85 

6.23 

5.7 

30.4 

NGC 2298 

4.75 

0.5 

5.88 

1.9 

15.3 

NGC 3201 

5.21 

-0.77 

5.90 

9.0 

22.1 

NGC 5053 

4.94 

-1.46 

5.38 



NGC 5139 

6.35 

-1.2 

6.91 

1.2 

6.2 

NGC 5272 

5.78 

-1.31 

6.30 

5.5 

13.4 

NGC 5286 

5.73 

-0.32 

6.58 



NGC 5466 

5.02 

-1.15 

5.52 

6.6 

57.1 

NGC 5904 

5.76 

-1.15 

6.33 

2.5 

35.4 

NGC 5927 

5.36 

-1.44 

5.81 

4.5 

5.5 

NGC 6093 

5.52 

-1.36 

6.02 

0.6 

3.5 

NGC 6121 

5.11 

-1.00 

5.75 

0.6 

5.9 

NGC 6205 

5.65 

-0.98 

6.29 

5.0 

21.5 

NGC 6218 

5.16 

0.1 

6.15 

2.6 

5.3 

NGC 6254 

5.22 

-1.1 

5.82 

3.4 

4.9 

NGC 6304 

5.16 

-1.85 

5.61 



NGC 6341 

5.52 

-1.23 

6.06 

1.4 

9.9 

NGC 6352 

4.82 

-0.6 

5.59 



NGC 6362 

5.01 

-0.49 

5.82 

2.4 

5.5 

NGC 6397 

4.88 

-1.6 

5.24 

2.34 

6.0 

NGC 6496 

5.12 

-0.7 

5.86 



NGC 6541 

5.64 

-1.07 

6.25 



NGC 6656 

5.63 

-1.5 

6.05 

2.76 

8.76 

NGC 6712 

5.23 

0.9 

6.53 

0.9 

6.2 

NGC 6752 

5.32 

-1.7 

5.59 



NGC 6809 

5.25 

-1.3 

5.77 

1.9 

5.8 

NGC 6838 

4.47 

0.2 

5.49 

4.8 

6.7 

NGC 7078 

5.90 

-1.9 

6.32 

5.4 

10.3 

NGC 7099 

5.21 

-0.92 

5.83 

3.0 

6.9 


erage a factor of 2 larger than theoretical predictions for 
the dynamical mass loss experienced b y Galactic clusters 
on their observed present- day orbits dKruiissen fc Mieskel 
2009; Rossi & Hurlcv 2015). Hence there is a discrepancy 
between our empirical calculation of the fraction of mass 
lost by a globular cluster using the present-day mass func¬ 
tion slope compared to the predicted fraction of mass lost 
from N-body simulations (for a given orbit). However it 
should be noted th a t we estimate initial masses less than 
IKruiissen fe Mieske] d2009h for NGC 288, NGC 5139, NGC 
6093, NGC 6121, NGC 6809, and NGC 7089. The latter two 
clusters having initial ma sses within 10% of the estimated 
IKruiissen fe Mieskel |2009l l values. This discrepancy, and the 
possible explanations behind it, are likely either due to real 
Galactic clusters not being subject to a comparable tidal 
field as our N-body simulations, or due to having initial pa¬ 
rameters outside the parameter space covered in this study 
(e.g. small initial c luster sizes or short perigalact ic distances 
). As discussed in lMarks, Kroupa fe Baumgardtl (|2008l ). pri¬ 
mordial gas expulsion may also help to unbind a significant 
fraction of low-mass stars in initially mass segregated clus¬ 
ters very early on in their evolution, and this effect is not 
accounted for in our simulations. 
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We know that our model clusters have not experienced 
the same tidal field as actual globular clusters since A) we 
have not modelled clusters evolving in their early formation 
environment and B) many clusters h av e orbits outsi de of our 
parameter space llDinescu et a l. 1999; Casct tFDinescu et al.l 
120071 . l2013l l . Including the additional mass loss mechanisms 
associated with a cluster’s early formation environment, 
specifically tidal shocks from encounters with GMCs and 
the unbinding of stars via gas expulsion, could explain why 
these clusters appear to have lost a significant amount of 
mass. However, we caution that the evolution of a prior to 
a cluster’s migration to the halo has not been explored, and 
may increase the uncertainty of the mass function - initial 
mass relation presented here. 

Modelling additional clusters with very small values of 
R p may partially explain the high initial masses of these 
clusters, as they are subject to an extremely strong tidal 
field and higher mass loss rates than the clusters m o delled 
in the present study. However iKruiissen fe Mieskel (120091 ) 
estimates that the mass loss rate due to tidal evaporation 
and disk shocking associated with clusters with small R p 
values still yield initial masses that are less tha n presented 
here. S imila r res ul ts were obta i ned b y Baumgardt fe Making 
(120031) and lBaumgardt, De Marchi fe Kroupal ( 20081 ). While 
a more detailed treatment of the Galactic potential at small 
R gc may partially decrease the discrepancy, it seems that 
the additional mass loss mechanisms discussed above will 
still be necessary to explain the high initial cluster masses 
associated with the empirically-measured stellar mass func¬ 
tions. 

Finally, it is possible that many observed clusters 
formed extremely compact compared to any of our model 
clusters, or were just recently placed on their current orbits. 
An extremely compact cluster will have time to almost com¬ 
pletely mass segregate before becoming tidally filling and 
would therefore have a higher present day a than our suite of 
simulations would predict. Similarly, if these clusters spent 
the majority of their lifetimes in a weak tidal field and were 
able to expand and relax before orbiting within the inner 
regions of the Milky Way, we would also observe a higher 
than expected a. 


5 SUMMARY 

We use a large suite of A-body simulations of star cluster 
evolution to determine the effects of initial mass, initial size 
and orbit on the evolution of a cluster’s stellar mass function. 
We illustrate that the observed slope of a cluster’s present 
day stellar mass function a is primarily dependent on the 
fraction of stellar mass that a cluster has lost and is approx¬ 
imately independent of its initial mass, initial size, orbital 
distance, eccentricity, and inclination. Our work confirms 
the results of previous numerical and analytical studies on 
the evolution of star cluste r s with long disrup tion time-scales 
llVesperini fe Heggiel 1 19971 : lKrmisseiJl2009bh . We note that 
there is some scatter about this relation due to the cluster or¬ 
bit, as clusters that are subject to stronger tidal fields, tidal 
heating or tidal shocks can be more easily s tripped of high 
mass stars, as found by IWebb et alj (l2014al) . There is also 
a secondary dependence of the stellar mass function’s evo¬ 
lution on the initial size of a cluster, since initially compact 


clusters have a chance to relax and mass segregate as they 
expand to fill their tidal radius. Due to the cluster under¬ 
going some degree of mass segregation before it starts to be 
tidally stripped, more low mass stars will be removed from 
the cluster and the slope of the mass function will increase 
faster. It also appears that orbital inclination can alter the 
evolution of a, as increased tidal heating due to passages 
through the Galactic disk can remove a larger fraction of 
higher mass stars from a cluster than if it orbited in the 
plane of the disk. 

Since the scatter induced by these factors is minimal, 
our models can be used to establish an analytic relation¬ 
ship between a and the fraction of mass that a cluster has 
lost (Equation [2]) that is accurate to within ~ 7%. If the 
initial size of a cluster can be constrained via other means, 
the uncertainties on the initial masses can be decreased to 
< 5%. However, we note that our mass function - initial 
mass relation is base d on the assumptions that all clusters 
form with a Universal lKroupa. Tout, fe Gilmor3 (1 19931) IMF 
and that forms of mass loss not included in our models, 
mainly tidal shocks experienced by clusters in their forma¬ 
tion environment, do not significantly alter the evolution of 
a. Additional simulations find that as long as the liigh-mass 
end of the IMF has —2.6 ^ ai ^ —2.3, our fitting func¬ 
tions are accurate. Clusters with top heavy IMFs tend to 
dissolve quickly with the stellar mass function undergoing 
very little change for stars between 0.15 — 0.8Mg. Clusters 
with bottom heavy IMFs lose very little mass due to stellar 
evolution such that a begins increasing much earlier than 
our mass function - initial mass relation predicts. Compar¬ 
isons to previous studies also suggest the maximum stellar 
mass and initial binary fraction may also influence the evo¬ 
lution of a. Future simulations are necessary to accurately 
model how these factors, as well tidal shocks experienced by 
clusters in their formation environment, will affect the early 
evolution of a. 

We use our mass function - initial mass relationship to 
determine the initial total stellar mass of any Galactic glob¬ 
ular cluster that has a known mass function and present day 
mass. Combining our mass function - initial mass r elation 
with the 33 Galactic globular cluster dataset of lPaust et al.l 
(l20ld) . we find surviving globular clusters were on average 
a factor of 4.5 times more massive than they are seen today, 
with some clusters being greater than a factor of 10 times 
more massive. We note that for two clusters our relation 
could only be used to put upper limits on their initial masses. 
We caution that these initial mass estimates are higher than 
other predictions in the literature based on mass loss due 
to tidal evaporation and disk shocking alone. This incon¬ 
sistency is also apparent from the A-body simulations we 
use to calibrate the relation between a and M/Mi (Equation 
[2J. If we extrapolate these simulations to globular cluster 
masses, they yield a total dynamical mass loss smaller than 
Equation [2] would predict for the observed mass function 
slopes of Galactic globular clusters. However, when scaled 
to globular cluster masses, our models would fail to repro¬ 
duce the observed distribution of stellar mass functions as 
they will have only lost a small percentage of their initial 
mass. This discrepancy suggests that additional mass loss 
mechanisms, like tidal shocks while in the early formation 
environment, are necessary to explain the present day mass 
functions of some Galactic globular clusters. Our predicted 
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initial masses would then only be in disagreement with pre¬ 
vious studies if these additional mass loss mechanisms do 
not alter the evolution of the mass function such that our 
predictions correspond to the initial mass of clusters once 
they have migrated into the halo. However if these additional 
mass loss mechanisms affect the evolution of a cluster’s mass 
function the same way that mass loss via tidally-limited two- 
body relaxation does, then our predicted initial masses ac¬ 
tually correspond to the true initial mass of each cluster at 
formation. The possibility still remains, however, that ad¬ 
ditional mass loss mechanisms may alter the evolution of 
the mass function in some non-standard way, in which case 
the connection between our predicted initial masses and a 
cluster’s true initial mass becomes unclear. The effects of 
additional mass loss mechanisms on the evolution of a need 
to be explored in order to determine how they influence the 
mass function - initial mass relation. 

Having a simple method to calculate initial cluster 
masses that is based on two observable parameters will help 
future studies to further constrain the GCIMF and initial 
distribution of cluster sizes. These are two very important 
parameters when studying both globular cluster and galaxy 
formation models. Our method provides a simple prescrip¬ 
tion for the evolution of a that can be adopted by fast codes 
wishing to model the evolution of large cluster popula tions 

dAlexander & Gieledl2012l . l2013l : [AIexander et al.ll2oH . e.g.) 

or be used to suggest initial parameters for simulations at¬ 
tempting to model individual globular clusters. 
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